Prom Exact to Partial Dynamical Symmetries: 
Lessons Prom the Interacting Boson Model 
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We exploit the rich algebraic structure of the interacting boson model to explain the notion of 
partial dynamical symmetry (PDS), and present a procedure for constructing Hamiltonians with this 
property. We demonstrate the relevance of PDS to various topics in nuclear spectroscopy, including 
K-bwad splitting, odd-even staggering in the 7-band and anharmonicity of excited vibrational bands. 
Special emphasis in this construction is paid to the role of higher-order terms. 
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INTRODUCTION 

The concept of dynamical symmetry (DS) is now widely accepted to be of central importance in 
our understanding of many-body systems. In particular, it had a major impact on developments 
in nuclear molecular [4] and hadronic physics [3], pioneered by F. Iachello and his colleagues. 
Its basic paradigm is to write the Hamiltonian of the system under consideration in terms of 
Casimir operators of a set of nested algebras. Its hallmarks are (i) solvability of the complete 
spectrum, (ii) existence of exact quantum numbers for all eigenstates, and (iii) pre-determined 
symmetry-based structure of the eigenfunctions, independent of the Hamiltonian's parameters. 

The merits of a DS are self-evident. However, in most applications to realistic systems, the 
predictions of an exact DS are rarely fulfilled and one is compelled to break it. More often one 
finds that the assumed symmetry is not obeyed uniformly, i.e., is fulfilled by only some states but 
C*"** ' not by others. The need to address such situations has led to the introduction of partial dynamical 

1 symmetries (PDSs) Q] . The essential idea is to relax the stringent conditions of complete solvability 

so that the properties (i)-(iii) are only partially satisfied. 

An exact symmetry occurs when the Hamiltonian of the system commutes with all the gener- 
^ ' ators of the symmetry group G. In this case, all states have good symmetry and are labeled by 

the irreducible representations (irreps) of G. The Hamiltonian admits a block structure so that 
incquivalent irreps do not mix and all eigenstates in the same irrep are degenerate. In making 
the transition from an exact to a dynamical symmetry, states which are degenerate in the former 
scheme are split but not mixed in the latter, in accord with the reduction G D G' D G" D 
and the block structure of the Hamiltonian is retained. Proceeding further to partial dynamical 
symmetry, some blocks or selected states in a block remain pure, while other states mix and lose 
the symmetry character. The hierarchy of broken symmetries and the corresponding Hamiltonian 
matrices are depicted schematically in Fig. 1. 
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FIG. 1: The Hamiltonian structure in an exact, dynamical (DS) and partial dynamical symmetry (PDS). 
In an exact symmetry, all states in a given irrep, on, are degenerate. The DS exhibits splitting but no 
mixing. In a PDS, only selected states (marked by an oval) remain solvable with good symmetry. 



IDENTIFYING HAMILTONIANS WITH PDS 

The algorithm for constructing Hamiltonians with PDS 0, Q starts from the chain of nested 
algebras 

Gdyn 3 G D ■ ■ ■ D G sym 
II I , (1) 

[h] (E) A 

where, below each algebra, its associated labels of irreps are given. Eq. (Q} implies that Gd yn 
is the dynamical (spectrum generating) algebra of the system such that operators of all physical 
observables can be written in terms of its generators; a single irrep of Gdyn contains all states of 
relevance in the problem. In contrast, G sym is the symmetry algebra and a single of its irreps 
contains states that are degenerate in energy. For N identical particles the representation [h] of 
the dynamical algebra Gd yn is either symmetric [N] (bosons) or antisymmetric [1 N ] (fermions) 
and will be denoted, in both cases, as [/ijv]. The occurrence of a DS of the type ([T]) signifies that 
the Hamiltonian is written in terms of the Casimir operators of the algebras in the chain, 7?ds = 
^2qClgC{G), and the eigenstates can be labeled as |[/ijv](E) ...A); additional labels (indicated 
by . . . ) are suppressed in the following. The eigenvalues of the Casimir operators in these basis 
states determine the eigenenergies £ , ds([^-at](E)A) of Hbs- Likewise, operators can be classified 
according to their tensor character under JT]) as Ti^mx- 

Of specific interest in the construction of a PDS associated with the reduction ([1]), are the 
rt-particle annihilation operators T which satisfy the property 

%„](a>A|[Mv£o)A)=0, (2) 

for all possible values of A contained in a given irrep (Eo) of G. Any n-body, number-conserving 
normal-ordered interaction written in terms of these annihilation operators and their Hcrmitian 
conjugates (which transform as the corresponding conjugate irreps) can be added to the Hamilto- 
nian with a DS (P), Hpns = #ds + J2 a A a p T^Tp, while still preserving the solvability of states 
with (S) = (Eq). The annihilation condition ^ is satisfied if none of the G irreps (E) contained 
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in the Gd yn irrep [/i/v-n] belongs to the G Kronecker product (a) x (So)- So the problem of 
finding interactions that preserve solvability for part of the states ([I} is reduced to carrying out a 
Kronecker product. 

In what follows we illustrate the above procedure and demonstrate the relevance of the PDS no- 
tion to nuclear spectroscopy. For that purpose, we employ the interacting boson model (IBM) [l], 
widely used in the description of low-lying collective states in nuclei in terms of N interacting 
monopole (s) and quadrupole (d) bosons representing valence nucleon pairs. The dynamical al- 
gebra is Gdyn = U(6) and the symmetry algebra is G sym = 0(3). The Hamiltonian commutes 
with the total number operator of s- and d- bosons, N = fi s + fid, which is the linear Casimir 
of U(6). The model accommodates three DS limits with leading subalgebras U(5), SU(3), and 
0(6), corresponding to typical collective spectra observed in nuclei, vibrational, rotational, and 
7-unstable, respectively. 



SU(3)-PDS AND AXIALLY DEFORMED NUCLEI 

The SU(3) DS chain of the IBM and related quantum numbers are given by [l| 

U(6) D SU(3) D 0(3) 
111. (3) 
[N] (X,(i) K L 

The multiplicity label K is needed for complete classification and corresponds geometrically to the 
projection of the angular momentum on the symmetry axis. The eigenstates \[N](X, fx)K, L) are 
obtained with a Hamiltonian with SU(3) DS which, for one- and two-body interactions, can be 
transcribed in the form 



Hds — h 2 



-G 2 (SU(3))+2N(2N + 3) + GG 2 (0(3)). (4) 



Here G 2 (G) denotes the quadratic Casimir operator of G. The spectrum of Hus is completely 
solvable with eigenenergies 

£ds = h 2 [-(A 2 +fi 2 + Xft + 3A + 3/i) + 2N(2N + 3)] + CL(L + 1) . (5) 

The spectrum resembles that of an axially-deformed rotor and the corresponding eigenstates are 
arranged in SU(3) multiplets. In a given SU(3) irrep (A,/i), each K- value is associated with a 
rotational band and states with the same L, in different i^-bands, are degenerate. In particular, 
the lowest SU(3) irrep (2N,0), contains the ground band g(K = 0) and the irrep (2N — 4,2) 
contains degenerate f3(K = 0) and j(K = 2) bands. This if-band degeneracy is a characteristic 
feature of the SU(3) limit of the IBM which, however, is not commonly observed in deformed 
nuclei. In the IBM framework, with at most two-body interactions, one is therefore compelled to 
break the SU(3) DS. 

To secure solvability and good SU(3) symmetry for the ground band, we follow the PDS algorithm 
and identify n-boson operators which annihilate all L-states in the SU(3) irrep (2N, 0). For that 
purpose, we consider the following two-boson SU(3) tensors, B^ n ^ x ^ K . em , with n = 2, (A, /i) = 
(0, 2) and angular momentum £ = 0,2 

4](0,2)0;00 « P = ^ ' ^ - 2 ( fit ) 2 > B [2](0,2)0; 2m « P L = ^ + V? ((ft rft)£> . (6) 
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FIG. 2: Spectra of 168 Er. Experimental energies (EXP) are compared with IBM calculations in an exact 
SU(3) dynamical symmetry [SU(3)], in a broken SU(3) symmetry (WCD) and in a partial dynamical SU(3) 
symmetry (PDS). The latter employs the Hamiltonian of Eq. ([9]), with TV = 16 and /12 = 4, 770 = 4, C = 13 
keV Q. 



The corresponding annihilation operators, Pq and p2 m , satisfy 

P \[N](2N,0)K = 0,L) = , 
P 2m |[7V](2iV,0)^ = 0,£) = , L = 0,2,4,...,2iV. (7) 

In addition, P satisfies 

P \[N](2N- Ak,2k)K = 2k, L) = , L = K, K + 1, . . . , {2N - 2k) . (8) 

For k > the indicated L-states span only part of the SU(3) irreps (A, n) = (2N — Ak, 2k) and form 
the rotational members of excited ^ k {K — 2k) bands. The combination PqPq+pJ '-^2 is completely 
solvable in SU(3) and is simply the first (/12) term in Hds, Eq. (@|, related to C*2(SU(3)). The 
SU(3)-PDS Hamiltonian is thus given by 

Ppds = Pds + m P^Po ■ (9) 

The relations in Eqs. ensure that Ppds retains solvable ground g(K — 0) and ^ k (K = 

2k) bands with good SU(3) symmetry (2N — Ak,2k) and energies as in Eq. ([5]). The remaining 
eigenstates of Ppds, including the (3{K = 0) band, are mixed. 

The SU(3)-PDS spectrum of the Hamiltonian © is compared with the empirical spectrum of 
168 Er in Fig. [5J As shown, the undesired /3-7 degeneracy is lifted and the PDS fit is of comparable 
quality to that of a broken-SU(3) calculation. Since the wave functions of the solvable states 0-© 
are known, one can obtain analytic expressions for matrix elements of observables between them. 
In particular, the resulting B(E2) values are found to be in excellent agreement with experiment Q, 
thus confirming the relevance of SU(3)-PDS to the spectroscopy of 168 Er. 

The dynamical symmetry expression, Eq. ([5]), implies a pure rotor spectrum with characteristic 
L(L + 1) in-band splitting. Such a pattern is observed in the empirical spectrum of the ground 
and (3 bands in 156 Gd, but the 7 band exhibits considerable odd-even staggering (OES). As shown 
in Fig. 3, the empirical staggering index, Y(L) = ^ 2L i ~ 1 ^ r^m^wjEgrj ~ 1, displays a pronounced 
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FIG. 3: Experimental (EXP) odd-even staggering in the 7-band of 156 Gd, compared with SU(3)-DS and 
SU(3)-PDS calculations. The latter employs the Hamiltonian of Eq. (JT5), with N = 12 and h 2 = 7.6, C = 
10.46, 772 = -0.24, 773 = 1.68 keV 0. The staggering index Y(L) is defined in the text. 



zigzag pattern, in marked deviation from a pure rotor for which Y(L) = 0. The fact that the SU(3)- 
DS is obeyed only in selected bands, highlights its partial nature. Following the PDS algorithm, 
we look for n-boson operators which annihilate the states in the g(K = 0) and (3(K — 0) bands. 
For n = 3, we identify the following SU(3) tensors, B},,^ ^ K . lm , with (A, /i) = (2, 2) and £ = 2, 3 

£[3](2,2)2;^ * = ^ = 2, 3 (10) 

which satisfy 

W em \{N}(2N,0)K = 0,L) =0 L = 0,2,4, ... , 2iV 

W £m |[7V](2iV-4,2)K = 0,i) =0 L = 0,2,4,...,(2jV-4) . (11) 

The Hamiltonian with SU(3)-PDS, now involving three-body terms, reads Q 

-ffpDS = ^ds + m Wl ■ W 2 + 773 Wl ■ W> 3 . (12) 

The cubic Casimir operator of SU(3), (?3(SU(3)), can also be included in i^DS- The relations in 
Eq. (fTTj) ensure that Hpns retains solvable ground g(K = 0) and (3{K = 0) bands with good SU(3) 
symmetry (2N, 0) and (2N — 4,2), respectively, and energies as in Eq. ([5]). Other eigenstates, 
including members of the 7 band, are mixed. As seen in Fig. 3, the SU(3)-PDS calculation can 
adequately reproduce the observed OES in 156 Gd [8(. In this case, the staggering arises from the 
coupling of the 7 band with higher excited bands. Other approaches advocating the coupling of 
the 7 band to the (3 band Q or to the ground band [l(]], cannot describe the OES in this nucleus. 
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0(6)-PDS AND 7-UNSTABLE NUCLEI 

The 0(6) DS chain of the IBM and related quantum numbers are given by [l| 

U(6) D 0(6) D 0(5) D 0(3) 
III I ■ (13) 

[N] (E) (r) n A L 

The multiplicity label n A is needed for complete classification. A completely solvable spectrum 
with eigenstates \[N]{T,)(r)n&L} and eigenenergies -Eds, is obtained with a Hamiltonian with 0(6) 
DS, which has the form 

H DS = h Q [-£ 2 (0(6))+#(# + 4)] + BC 2 (0(5))+CC 2 (0(3)) , 

Eds = fro [-E(E + 4) + N(N + 4)] + B t(t + 3) + C L(L + 1) . (14) 

The spectrum resembles that of a 7-unstable deformed rotor, where the states are arranged in 
bands with 0(6) quantum number E = N — 2v, (v = 0,1,2,...). The ground band (v = 0) 
corresponds to the 0(6) irrep with E = N. The 0(5) and 0(3) terms in 7?ds (O, govern the 
in-band rotational splitting. 

The 0(6)-DS limit provides a good description of the empirical spectrum and E2 rates in 196 Pt, 
for states in the ground band (E = N). This observation was the basis of the claim [11] that the 
0(6)-DS is manifested empirically in this nucleus. However, as shown in Fig. 2J the resulting fit to 
energies of excited bands is quite poor. The 0^, O3 , and O4 levels of 196 Pt at excitation energies 
0, 1403, 1823 keV, respectively, are identified as the bandhead states of the ground (v = 0), first- 
(v = 1) and second- (v = 2) excited vibrational bands [11]. Their empirical anharmonicity, defined 
by the ratio R = E(v = 2)/E(v = 1) - 2, is found to be R = -0.70. In the 0(6)-DS limit these 
bandhead states have r = L = and E = N, N — 2, N — 4, respectively. The anharmonicity 
R — —2/(N + 1), as calculated from Eq. (TT41 . is fixed by N. For N = 6, which is the appropriate 
boson number for 196 Pt, the 0(6)-DS value is R = -0.29, which is in marked disagreement with 
the empirical value. Large anharmonicities can be incorporated in the IBM only by the inclusion 
of at least cubic terms in the Hamiltonian [12J . One is therefore confronted with the need to select 
suitable higher-order terms that can break the DS in excited bands but preserve it in the ground 
band. These are precisely the defining properties of a PDS. Following the general algorithm, we 
look for n-boson operators which annihilate all the states in the 0(6) ground band. For n = 3, 
there are two such 0(6) tensors, -B^^^^.^, with a = 1, (t, i) = (0, 0) and (r, £) — (1, 2), given 
by 

^[3]<l)(0)0;00 X P Q S ^ ' ^[3](l)(l)0;2m a 44n ■ (15) 

Here Pq = & -S — (s^) 2 , is an 0(6) scalar and PqPo is simply the first (ho) term in -Hds, Eq. (fl4")) . 
related to (72(0(6)). The operators, sPo and d m Po, annihilate all (t, L)-states in the 0(6) irrep 
(E) = (N). This is ensured by the fact that 

P \[N]{N)(T)n A L) =0 , t = 0,1,2,...,JV. (16) 

The only three-body interactions that are partially solvable in 0(6) are thus Pjn s Po and Pgn^Po. 
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FIG. 4: Observed spectrum of 196 Pt compared with the calculated spectra of Hds (|14p . with 0(6) dy- 
namical symmetry (DS), and of Hpos (|17[) with 0(6) partial dynamical symmetry (PDS). The parameters 
in Hns (Hpds) are ho = 43.6 (30.7), B = 44.0 (44.0), C = 17.9 (17.9), and p = (8.7) keV. The boson 
number is N — 6 and E is an 0(6) label 



Since the combination PjC^s + )P = {N - 2)P f Po is completely solvable in 0(6), we can 
transcribe the 0(6)-PDS Hamiltonian in the form 

-ffpDS = ^ds + Po PfaPb ■ (17) 

The spectrum of Hpr>s (fT7|> is shown in Fig. |4j The states belonging to the S = N = 6 multiplet 
remain solvable with energies which obey the same DS expression, Eq. fH|) . States with £ < 6 
are generally admixed but agree better with the data than in the DS calculation. For example, 
the bandhead states of the first- (second-) excited bands have the 0(6) decomposition S = 4: 
76.5% (19.6%), E = 2: 16.1% (18.4%), and S = 0: 7.4% (62.0%). Thus, although the ground band 
is pure, the excited bands exhibit strong 0(6) breaking. The calculated 0(6)-PDS anharmonicity 
for these bands is R = —0.63, much closer to the empirical value, R = —0.70. It should be 
emphasized that not only the energies but also the wave functions of the S = N states remain 
unchanged when the Hamiltonian is generalized from DS to PDS. Consequently, the E2 rates for 
transitions among this class of states are the same in the DS and PDS calculations 0]. Thus, 
the additional three-body (po) term in the Hamiltonian (TT71) . does not spoil the good 0(6)-DS 
description for this segment of the spectrum. 



SUMMARY AND CONCLUSIONS 

The notion of partial dynamical symmetry (PDS) extends and complements the fundamental 
concepts of exact and dynamical symmetries. It addresses situations in which a prescribed sym- 
metry is neither exact nor completely broken. When found, such intermediate symmetry structure 
can provide analytic solutions and quantum numbers for a portion of the spectra, thus offering 
considerable insight into complex dynamics. 

On phenomenological grounds, having at hand a concrete algorithm for identifying and con- 
structing Hamiltonians with PDS, is a valuable asset. It provides selection criteria for the a priori 
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huge number of possible symmetry-breaking terms, accompanied by a rapid proliferation of free- 
parameters. This is particularly important in complicated environments when many degrees of 
freedom take part in the dynamics and upon inclusion of higher-order terms in the Hamiltonian. 
In the IBM examples considered, there are 17 possible three-body interactions, yet only a few 
terms satisfy the PDS requirements. Futhermore, Hamiltonians with PDS break the dynamical 
symmetry (DS) but retain selected solvable eigenstates with good symmetry. The advantage of 
using interactions with a PDS is that they can be introduced, in a controlled manner, without 
destroying results previously obtained with a DS for a segment of the spectrum. These virtues 
greatly enhance the scope of applications of algebraic modeling of quantum many-body systems. 
On a more fundamental level, PDS can offer a possible clue to the deep question of how simple 
features emerge from complicated dynamics. 

PDSs are not restricted to a specific model but can be applied to any quantal systems of interact- 
ing particles, bosons, as demonstrated in the present contribution, and fermions 13-15J. They are 
also relevant to quantum phase transitions [16| and to the study of mixed systems with coexisting 
regularity and chaos [3, [l 
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